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Abstract. In this article we analyze the radiation loss from a high energy cosmic 
ray proton propagating in a spacetime with non-systematic Lorentz violation. From 
an effective field theory perspective we illuminate flaws in previous attempts that 
l/^ I use threshold approaches to analyze this problem. We argue that in general such 

approaches are of rather limited use when dealing with non-systematic Lorentz 
violating scenarios. The main issues we raise are a) the limited applicability of 
ly^ ' threshold energy conservation rules when translation invariance is broken and b) 

■ the large amounts of proton particle production due to the time dependence of the 
fluctuations. Ignoring particle production, we derive a constraint on the magnitude of 

■ velocity fluctuation \vf\ < 10~^-^, much weaker than has been previously argued. 
lO ■ However, we show that in fact particle production makes any such constraint 

completely unreliable. 
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+i . 1. Introduction 

a 

Testing quantum gravity prima facie seems to be an impossible task, since the natural 
^ I scale of the theory, the Planck energy of 10^^ GeV, is sixteen orders of magnitude 
' in energy beyond the reach of earth-based accelerators. If we cannot somehow gain 
knowledge about Planck scale physics, however, then a complete theory of quantum 
gravity will most likely remain out of reach. Fortunately, even if it is not possible to 
directly probe Planck-scale physics, quantum gravity may leave tiny residual effects at 
lower energies that are observationally testable. A promising avenue in this vein, and 
certainly the one that has received the most attention, is the possibility that Lorentz 
invariance is violated at the Planck scale. 

If Lorentz symmetry is broken by quantum gravity, then it is only an approximate 
symmetry and there should be small Lorentz violating effects in low energy physics. 
There are two main approaches to modelling low energy Lorentz violating effects. The 
first is the standard model extension p]. In this approach one follows the tenets of 
effective field theory and adds all possible Lorentz violating, renormalizable operators 
to the standard model. This definite (and large) set of operators yields a number of 
different effects, such as sidereal signals in clock comparison experiments c.f. [2] or 
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modified neutrino oscillations jB], which can be used to constrain the Lorentz violating 
coefficients in the standard model extension. For a description of the current bounds 
on the standard model extension see |1] and references therein. 

An alternative approach is that of modified dispersion for elementary particles, 
since one would expect Lorentz violating corrections to the usual particle dispersion 
E"^ = + 777,2,11 Usually in the literature rotation invariance is assumed to hold in some 
frame, for the following reason. If the rotation subgroup of the Lorentz group is broken 
but the boost subgroup is preserved, then rotation breaking would occur at all energies. 
Since rotation symmetry is a very good symmetry at low energies, Planck scale rotation 
breaking is either tiny or accompanied by broken boost invariance. The same is not 
true of a violation of boost symmetry. We have explored only a tiny fraction of the 
boost subgroup, which is after all unbounded. Hence it is possible that boost Lorentz 
violation could be small at low energies but become significant as the energies approach 
the Planck scale. From a phenomenological standpoint is is therefore logical to look for 
broken boost invariance first while preserving rotation invariance. 

Under the assumption of rotation invariance a dispersion modification would 
typically look hke 

E'=p' + m\ + r^A^ (1) 

where rj^ is usually taken to be an 0(1) coefficient that can depend on particle species 
A. rjA can also depend on particle properties such as helicity, however we shall 
ignore such possibilities as they are not relevant for our analysis, n is an integer 
determined by the underlying assumptions about quantum gravity. Note that if n > 2 
then these dispersion modifications would come from non-renormalizable terms in a 
particle Lagrangian. Modified dispersion affects the kinematics of reactions, leading 
to anomalous thresholds — reactions usually forbidden can be allowed, and threshold 
energies for already allowed reactions can be raised or lowered EI • The usefulness of 
modified dispersion is that the energy at which the reaction kinematics change is often 
far below the Planck energy. 

To see this, consider the effect of Lorentz violation with n = 4 on the existence of 
ultra-high energy cosmic rays (UHECR), as done in jHlEllSI- We assume for this analysis 
that the UHECR are actually protons. If rjp > then protons above a threshold energy 
of approximately Eth ~ {m'^MpJripY/'^ are superluminal and lose energy rapidly via 
Cerenkov emission. In such a scenario there would be a sharp break in the UHECR 
spectrum at Eth as there could be no long lived protons above Eth- For an rjp of 0(1), 
Eth ~ 10^^ eV. Since we observe a UHECR spectrum out to 10^° eV, rjp is constrained 
to be much less than 0(1). So we see that we can constrain Lorentz violating terms 
with current astrophysical measurements even if those terms are Planck suppressed, as 
expected if they arise from quantum gravity. 

II An effective field tfieory of course yields particle dispersion relations. However, often a dispersion 
modification alone is studied without explicitly working in a particular field theory. 
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A modification of tlie form in eq. (1) is an example of systematic Lorentz violation, 
which is constant over the lifetime of a particle. However, some ideas about quantum 
gravity [3 El E] give rise to stochastic fluctuations in the world-line of a particle over 
its lifetime. The question then arises, what happens if we combine the two ideas of 
Lorentz violation and non-systematic effects? In this paper we are concerned with 
the phenomenology of such non-systematic Lorentz violation, especially with the non- 
systematic version of the Cerenkov effect.^ There has been some previous work on non- 
systematic scenarios and their possible observational signatures. Amelino-Camelia ^Hj 
has looked at possible imprints of metric stochasticity on gravity wave interferometry. 
The basic conjecture is that fluctuations of geometry induces a new noise-strain that is 
within the range of detectability of space-based observatories such as LISA. However, 
for this to happen it must be assumed in the model that fluctuations add up coherently 
along the entire length of the interferometer arms,"*" which is a very fine tuned scenario 
that seems unlikely. 

Aloisio et. al fTJ |E| and others fH^ ITi] have considered how non-systematic 
dispersion fluctuations affect the propagation of ultra-high energy cosmic rays (UHECR) 
and analyzed shifts in the thresholds for processes such as the GZK reaction, p + 7cmb — 
p + vr*^. We shall focus on the conclusions of Aloisio et. al., although the general 
arguments are applicable to any non-systematic threshold analysis. In ^TJ ^] the 
authors use a non-systematic dispersion relation similar to (P) where the coefficient r/A 
varies randomly as the particle propagates. It is further assumed that the timescale 
of the variation is set by the de Broglie frequency 1/E of the UHECR. With these 
assumptions they then derive that the GZK-threshold shifts to sub-GZK regimes and 
conclude that observations below the GZK energy scale are enough to constrain scenarios 
of fluctuating spacetimes. Furthermore they find a stability problem with UHECR, 
namely that all charged particles are unstable to the vacuum Cerenkov emission of 
photons. This conclusion, which is a primary focus of this paper, seems encouraging 
at first sight, since sub-GZK cosmic rays follow a simple power law spectrum and are 
obviously stable over cosmological times. Fluctuating models as considered by Aloisio 
et. al. would therefore seem to be ruled out. 

The analysis of Aloisio et. al. only considers the effects of fluctuating dispersion 
relations on energy conservation. There is a problem with this sort of analysis however - 
the length scale of dispersion fluctuations may be shorter than the de Broglie wavelength 
of other particles in the reaction (as pointed out in PS])- There are many fluctuations 
in the interaction region and hence no single energy-momentum conserving (5-function. 

Furthermore, fluctuations require some time-dependent background field (usually 
the metric) to be present. A quantum field in a time dependent background will 
generically experience particle production - a competing effect to radiative energy 
loss. In fact depending on the assumptions made, the energy loss via radiation can 
be completely overwhelmed by particle production of UHECR protons. The resulting 

% For a discussion of the systematic Cerenkov effect, see for example 
+ We thank Steve Car hp for pointing this out. 
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spectrum of UHECR in such scenarios would still have a flux at very high energies, 
even though each individual proton is radiating and evolving towards low energy, as 
the time dependent background pumps UHECR protons into the spectrum at very high 
energies. Constraints based on the absence of UHECR protons above a certain energy 
must therefore take this particle production into account. 

We have two primary goals. The first is to illuminate in more detail why simply 
deriving thresholds from an energy-momentum conservation equation and using those 
thresholds to establish constraints can be misleading. Secondly, we construct a toy 
field theory model and analyze the energy loss and particle production rates for high 
energy cosmic rays to estimate the magnitude of any constraints that can realistically be 
established. These constraints turn out to be quite weak (at best). Our focus will be on 
the non-systematic analogy of the vacuum Cerenkov effect, although some of the general 
conclusions will be applicable to any reaction involving low energy particles (such as the 
GZK reaction). In more detail, section |21 discusses the problems with the rate calculation 
of Aloisio. In section El we build a model for a representative UHECR proton (from a 
population at some energy) as it propagates in a stochastic background, calculate the 
radiation spectrum of this model, and explore the implications for UHECR's. We discuss 
particle production and its associated complications in section HJ Finally, in section El 
we summarize and speculate about possible future directions. Throughout this paper 
we adopt units such that h = c = 1 and choose metric signature + — — — . 



2. Non-systematic dispersion and energy thresholds 

We start by reviewing the approach taken by Aloisio to derive the charged particle 
instability from the vacuum Cerenkov effect p — > p + 7. In some preferred frame, 
usually taken to coincide with the rest frame of the CMBR, the dispersion relations 
assumed by Aloisio et. al. are 

E'=p' + m' + (2) 



E'' = + m' + (3) 
k' + V.^ (4) 



/3 



UJ 



.2 _ ; 2 



Mpi 



where p,p',k are the momenta of the incoming proton, outgoing proton, and photon 
respectively. The coefficients rip, rj^ are taken to be independent and randomly chosen 
according to a gaussian distribution with some amplitude over every de Broglie 
wavelength of the particle. The authors motivate this independence by noting that 
the time scales involved in the interaction are all much greater than the Planck time 
and therefore particles should feel different underlying quantum gravity fluctuations. 
These modified dispersion relations are then substituted into the conservation equation 

P/. - = P'f.- (5) 
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Furthermore, the authors argue that since the fluctuations are independent the outgoing 
particle fluctuations can be ignored without significantly changing the results (as for 
a significant fraction of the time the outgoing particles will have Lorentz violating 
coefficients much less than the incoming proton). Squaring both sides of (0) then leads 
to an equation for the threshold momentum where the reaction is kinematically allowed, 



In ® we have taken all the particle momenta to be parallel as this is the appropriate 
configuration for a threshold jT2j. As u —>■ 0, the threshold momentum also goes to 
zero, independent of the size of rjp, which implies that charged particles would always 
be radiating low energy photons. This is the instability as discussed in Aloisio et. al. 

We now examine the validity of this procedure. Implicit in this approach is the 
assumption that the reaction rate is fast once kinematically allowed so that the relevant 
quantity for a constraint is the threshold. In the systematic vacuum Cerenkov effect 
this is true, the energy loss rate can be shown to rapidly approach E^/Mpi for momenta 
only shghtly above pth- A GZK proton with energy 10^^ eV would then radiate most 
of its energy on the order of 10~^^ seconds. However, the same assumption cannot be 
made in the non-systematic case, especially for low energy emitted photons. In the 
usual approach to scattering in quantum field theory the 5-function constraint 
that enforces energy conservation is a consequence of the adiabatic approximation — 
namely that the interaction potential (the time dependent background geometry causing 
the fluctuations in this case) varies slowly over the region of interaction. With the 
Aloisio et. al. assumption that the fluctuation time scale t/ is set by the energy of 
the particle then ~ much shorter than the wavelength of an emitted low energy 
photon. This immediately implies that a low energy photon will interact with many 
fluctuations, thereby invalidating the use of a (5-function that enforces conservation 
with one particular fluctuation. More severely, however, the adiabatic approximation 
fails and one will not in general have a conserved energy. This invalidates the whole 
notion of a simple threshold analysis. Instead, there will be particle production which 
must be taken into account in addition to any radiative losses. 

If the random choice of the coefficients rip happens to add up coherently over the 
interaction region, then a kinematic threshold analysis is legitimate, as the interaction 
potential is a constant over the interaction region. The question then becomes, is this 
likely over the lifetime of a UHECR with emission of a low energy photon? Even if the 
threshold formalism does not apply in general there may be a high probability for it to 
apply at least once and hence we might still have an instability. We now estimate this 
probability in a simple model. First, let us consider the simplest possible fluctuation, 
where rjp takes values ±1 on a timescale of tf = 1/E. For a UHECR of energy 10^® 
eV that travels over 1 Gpc, there are 10^^ total fluctuations. If the fluctuations are to 
add up coherently for some timescale given by l/uj, then for u = E/n where n > 1 the 
probability is 2~". Applying some simple probability shows that for n > O(IOO) there is 
almost zero chance that the fluctuations will add up coherently at any point during the 




(6) 
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UHECR lifetime. This implies that a threshold approach might only apply if we restrict 
the outgoing photons to have energies above 10^^ eV. Note that even this number is 
conservative as we have forced the fluctuations to have only two possible values, thereby 
increasing the probability that they might be coherent. 

Even in the high energy photon region of phase space we have an issue, however. 
As Aloisio et. al point out, it is also assumed in their analysis that the fluctuations in 
each of the particles' energy is independent of that of the others in the spacetime region 
of interaction — that is each of the coefficients tja have a different independent value. As 
seen however, a threshold analysis is only applicable when dispersion fluctuations act 
coherently over the entire interaction region. If this is the case, it is quite reasonable that 
the underlying spacetime fluctuations that give rise to a coherent dispersion fluctuation 
might also be coherent. The effective length scale of the spacetime fluctuation would 
therefore be on the order of the interaction region. It seems suspect to treat each 
dispersion modification as independent in such as scenario, as the underlying spacetime 
fluctuation is the same for all particles. In short, we doubt that in a threshold 
analysis the Lorentz violating coefficients can be considered independent. Unfortunately, 
knowledge of how they are correlated would require a more specific model. If the 
dispersion fluctuations were independent, a threshold analysis similar to that of Aloisio 
might be valid, although it would necessitate restricting the outgoing photon to be high 
energy. Such an analysis has not been performed; we shall not do so here since we doubt 
an independent coefficient scenario is applicable. 

So what can one do? Note that the discussion above did not argue that field theory 
fails, rather merely the implementation in terms of a single energy-conserving 5-function 
in the rate is not applicable. We can still analyze the emitted radiation, we just need to 
use the appropriate field theory techniques. In the next section we estimate the energy 
loss rate for low frequency radiation (to < E/lOO), assuming that the UHECR dispersion 
fluctuations are incoherent. This will allow us to show that there is in fact no stability 
problem for UHECR protons. We use classical field theory for this estimation since we 
are restricting to low energy outgoing photons (we can neglect any individual photon's 
backreaction on the UHECR). 

3. Radiative emission in a fluctuating background 
3.1. General Considerations 

The vector potential from a moving particle (proton) p with charge e is given in Lorentz 
invariant electromagnetism by the expression 



where Dr{x — x') is the usual retarded Green's function. J°'{x') is the current, given in 
covariant form by 




(7) 




(8) 



Constraints from cosmic rays on non- systematic Lorentz violation 



7 



where r is the proper time,f"(r) is the four- velocity and r(r) is the proton's position. 
There are two ways that fluctuations can enter into (|7j): the Green's function or the 
current. We first deal with the Green's function. 

One would naturally expect that spacetime fluctuations would modify the Green's 
function for the electromagnetic potential A°' (the equivalent of a the modification of 
the photon dispersion in the previous section). Unfortunately, by definition in Lorentz 
violating field theories there are couplings not only to the metric, but also to the 
fields that are responsible for Lorentz violation. For example, in the aether model 
of pni; boost invariance is broken by the introduction of a new unit time-like vector 
field u°' which can couple to via terms such as u^u^ FapF^ , Mpl F^a{u^ df3)u'^ , 
or Mpfu^u^Fai3{u^ds)'^F^. New couplings in combination with the fact that the 
fluctuations yield time dependent values for the metric and u°' or other Lorentz violating 
fields make it nigh impossible to evaluate the electromagnetic potential with a complete 
Lorentz violating, fluctuating Green's function. We can still make progress if we assume 
that Lorentz violation scales with energy, as in the case of the dispersion relations in 

If this is the case, then since we are looking at low frequency emission relative 
to the UHECR energy, any Lorentz violating fluctuations in the electromagnetic sector 
will be much smaller than the fluctuations in the source current. Hence, for the rest of 
this discussion we will ignore Lorentz violation in the electromagnetic sector and only 
consider the usual Lorentz invariant Green's function for the electromagnetic field. 

We now turn to the question of Lorentz violation in the current J". Since the 
current is only a function of position (directly and through the velocity) there is actually 
no need to use the fluctuating dispersion framework that has been previously preferred. 
Instead, we will consider the effect of a fluctuating velocity and relate it to dispersion 
relations only when comparing with other constraints. The largest radiation response 
occurs for the component of the velocity fluctuations in the direction of propagation 
(which we choose to be along the z-axis), so we will look only at this case. Finally, 
we must choose some model for the fluctuations. We shall choose random velocity 
fluctuations in the z-direction with values ±Vf and some timescale t/ as this can be 
directly related to the fluctuating dispersion scenario. More sophisticated choices are 
of course possible, for example vj could be chosen from some underlying probability 
distribution. We choose the simple bivalued random walk scenario as it captures the 
essential physics without added complications. In summary, we choose our charged 
particles to have a velocity in the z-direction with magnitude v = vq ± Vf, where the 
sign is randomly chosen every tf. 

The deviation from the Lorentz invariant trajectory of any particle in a population 
will be different, obviously, if each particle experiences independent fluctuations. This 
is extremely likely, since our observed population of UHECR's has sources that vary 
over cosmological distances. From an observational standpoint it is irrelevant if a single 
UHECR proton out of a large population experiences a significant amount of energy 
loss from radiation. The only possible observational signature is if a large fraction of 
protons lose energy from fluctuations, leading to a cutoff in the spectrum above some 
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energy. A population with some initial Vq (alternatively at some energy E) will spread 
out over time in z around the value ZQ{t) = v^t. Since we have chosen the random walk 
scenario for fluctuations, the variance of the population is simple to calculate. Each 
particle deviates in the z-direction hy Zf = ^Vftj every step and in a time interval t 
there are t/tf steps, which yields a variance in z of Vf ^Jlt^tj-K. We therefore use as a 
representative trajectory 



z{t)={v^t^VfS^-t^t) (9) 

i.e. a particle with a deviation from the Lorentz invariant trajectory equal to the variance 
of the population (assuming it was created at t = 0). If a particle with this trajectory 
radiates much of its energy, then we know that a large fraction of the population is 
radiating and vice versa. Similarly, the velocity corresponding to Q is 



Equations and (jiup determine the current J". Note that some of the particles will 
slow down while others speed up since the sign in front of can be either positive or 
negative. 

As an aside, we mention that since z{^\v completely determine the representative 
current, our results can be quickly be generalized to other distributions. Models with 
anti-correlations between successive fluctuations have representative trajectories that 
deviate much more slowly than the random walk model we are considering. In general, 
if the rate of deviation of z{t) versus time grows more slowly than v^, the corresponding 
energy loss will be less (and vice versa if the deviation grows more quickly). Hence the 
constraints from energy loss on these models are presumably weaker. In particular the 
"holographic" model of Ng and van Dam [H] would be less constrained. For a discussion 
of various models for the deviation of particle trajectories from spacetime foam see |2Uj . 

3.2. Energy loss rate 

Since we are considering the effect of Lorentz violation only in the current J", we can 
easily calculate the energy loss rate. The energy radiated per unit frequency and solid 
angle is |2I] 

d^l _ 
dujdVL 47r2 

where n is the unit normal between x and the observation point. Given and (jlOj) 
the 3-current for a particle with charge e is 



dt d^xnx[nx J(x,t)]e^"(*~"-^") 



2 

(111 



J{x,t) = ze(i;o±i;/y^)5(x)5(y)5(z-t;otT^/y|t^)e(t)e(T-t).(12) 

The Heaviside step functions Q(t),Q(T — t) are present as we are assuming the particle 
is created by some astrophysical source at time t = and observed on earth at time 
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t = T. Substituting (jl2j) into (jllj) and working in spherical coordinates we have 

-sin'^9 



dt (t;o±t;^W^)e*-*(i-"o™^^^''/=-^\/^) (13) 



where A = 1 — i;ocos^^. The integral in p3|) can be evaluated explicitly, 



dt (t;o ±t;j-^/^)e^"*(^-^°'°'^^''/™'^V^^) 



2e 4A 
\AJA 



R 



( sin(x^) — cos(x^)) 



x=(3 



x=7 



(14) 



where we have defined 

B = ±Vf cos 9 

13 = VujAT - 
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R 



UJ 
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4A 



B 



IujA 



(15) 



C and S are the Fresnel cosine and sine functions respectively, defined by C{u)\iS{u) 



JJ^ga'^y (iy. Combining all these elements, the radiated energy is given by 



d^I 



ijjt 



sin^^ 



dwd'^ n'^A 
+— ( sin(x^) — i cos(x^)) 



(16) 



x'=/3 



x=7 



3.3. Analysis 

The radiated energy (fT^ is not all due to the effect of the fiuctuations. If we set f / = 0, 
then there is a residual piece 

^ = ^--^--^ (17) 

which is due to the finite time existence of the source. We will neglect this piece since 
we are looking for the energy radiated during the particle's flight, not the energy needed 
to create or destroy the source particle. The other terms are fluctuation dependent and 
reflect radiation loss over the entire travel time. 

The remaining frequency and angular integrals in (jl6|) are extremely complicated 
analytically but numerically tractable (although difficult). The numeric difficulty arises 
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because /3,7 can be large, which makes the Fresnel functions and the sin(x^), cos(x^) 
terms rapidly oscillatory. We deal with this situation by analytically replacing the 
appropriate Fresnel function with its asymptotic limit in the relevant regions of phase 
space. We then rewrite the integrand averaging over the oscillations, i.e. we replace 
sin^ , sin^ 7^ terms by 1 /2 and set cross terms such as sin 7^ sin to zero as they are 
almost orthogonal. This will introduce small errors, however this approximation will 
not significantly affect our results. 

It is only possible to give the total energy loss as a function of A/ = ±Vfy/tf, as 
this combination of parameters completely controls the energy loss in (fT^ . Hence we 
can only constrain Af and not Vf,tf individually. Since we have assumed that Lorentz 
violation must scale with energy, the observationally most sensitive case is at the high 
end of the UHECR spectrum, the near GZK protons, as they have both the largest 
known energy and gamma factor. Furthermore, they are presumed to be extragalactic 
and so travel over cosmological distances, enhancing the effect of any small energy loss. 
Hence for our energy loss estimate we take a proton of energy 10^^ eV travelling over a 
distance of 1 Gpc. The energy loss as a function of is a power law, best fit by the 
curve 



There is a small deviation of a fraction of a percent in the energy loss between the cases 
Af > and Af < 0. The difference is low because the energy loss for our parameter 
choice is dominated by the term in (fT^ which goes as (±vj)^. The important 
physical effect is that protons that slow down or speed up emit in almost the same way. 
We will neglect the small difference when deriving constraints. 

If we temporarily assume, as has been done previously, that there is no particle 
production of UHECR's from the fluctuations then we can derive constraints on Vf,tf 
in the following manner. We observe UHECR protons at energies of 10^^ eV (and above). 
If these protons are created at the source with an energy near 10^^ eV, then the total 
energy loss during their flight must be much less than 10^^ eV or else we would not see 
any flux in that energy band. Demanding that E^ad < 10^^ eV yields the constraint 
\Af \ < 10~^^'^. It is possible, although unlikely, that high energy protons are created 
with higher energies, lose energy during flight, and reach earth at the observed energies. 
If we are extremely conservative and take the worst case scenario, that the protons 
that reach us at GZK energies are created at the Planck energy, then the energy loss 
can at most be the Planck energy. This yields the constraint |A/| < 10^^^'^. We shall 
not consider this possibility, instead choosing the much more plausible scenario where 
UHECR cosmic rays have energies at the source near their observed energies. 

The constraint |A/| < 10~^^-^ can be used to put constraints on Vf, assuming some 
length scale for tf. There are two length scales intrinsic to the problem, the Planck 
length and the de Broglie wavelength of the UHECR proton. If tf = 1/ E, a.s in Aloisio 
et. al. then Vf < 10~^'^. This constraint is very weak in comparison with other results, 
as we show in section EiH If tf = tpianck then the constraint is further weakened to 




(18) 
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Vf < 10~^'^, however we remind the reader that we have derived these constraints in 
a classical framework. If the fluctuation time is less than the de Broglie wavelength 
of the UHECR proton presumably the quantum nature of the UHECR must be taken 
into account. We include the tj = tpi case simply to show that the constraint weakens 
as tf — > tpi. This failure does not affect our main conclusion, that there is no way to 
realistically calculate useful constraints using this method. 

The constraints above only take into account frequencies up to E/lOO = 10^^ 
eV. For informative purposes we show below the corresponding expressions allowing 
fc = 10^^, 10^^ eV. These expressions should not be considered constraints, since they 
might suffer from the correlation problems discussed above. However, they show that 
even if we push all the way to 10^^ eV, the constraints would be improved by not 
even two orders of magnitude. 

fc = lO^^eV : log Erad = 2.07% |A/| + 54 (19) 
fc = lO^^eV : log Erad = 2.06% |A/| + 54.9. (20) 

(21) 

3.4- Comparison with other constraints 

We now compare the strength of our constraint, vj < 10^^'^, with other direct 
constraints on particle velocity fluctuations. We first contrast our results with those 
for systematic Lorentz violation. The existence of stable UHECR protons at energies of 
10^^ eV requires that the speed of protons does not exceed the speed of light by more 
than one part in 10^°, i.e. 

Vp-v^< 10-20 (22) 

which is thirteen orders of magnitude stronger than our limit onvf. Note that our limit 
is, however, two-sided as both positive and negative fluctuations emit radiation. 

If we assume that the velocity fluctuations come from a modified dispersion relation 
via V = dE/dp then the velocity fluctuations can be written as 

= " ^^^^ (23) 

where rjp is the 0(1) coefficient for a proton. Our corresponding constraint on r]p is 

2 lO"^'^ 

< ^_l(10-9)n-2- (24) 

For n > 2, as we must have since we assume Lorentz violation scales with energy, we 
see that the resulting constraints are far greater than 0(1), rjp < 0(10^), 0(10^^) for 
n = 3,4 respectively. Hence with fluctuating models there are no strong constraints that 
can be placed on n = 3,4 fluctuating dispersion, contrary to the conclusions in |llj. 

Furthermore, our results show that there is no particle stability crisis. Our 
calculation includes the contribution to the energy loss rate from low frequency radiation 
and shows that it is actually finite and small. It is true (as one can see from |TB])) that 
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given any fluctuation in f j for any energy proton tliere is some radiation. This rougfily 
corresponds to tlie statement in for any value of Pin,Vp there is some u that is at 
threshold, i.e. some radiation. However, the loss rate is small enough that there is no 
stability problem for a cosmic ray proton over a cosmological travel time. 

4. Particle production 

In the above analysis, we have restricted ourselves to fluctuations in the velocity of a 
particle, which is enough to calculate the radiated energy classically. Previous work has 
not taken into account the quantum nature of the UHECR protons, and in particular the 
consequences of particle production. Without particle production a continuous energy 
loss will actually deplete the population of high energy protons, which is a necessary 
process for the observational constraint. Unfortunately, in a fluctuating background 
this method does not work. Systematic Lorentz violation * preserves translation 
symmetry - there exists a global Killing vector by which a conserved particle energy and 
number can be defined. If the underlying background is fluctuating, however, no such 
Killing vector exists, time-translation symmetry is broken, and there will generically be 
particle production. One should expect significant particle production in the Aloisio 
et. al. model since the population of UHECR is being driven by the fluctuations on 
a time scale of order of their de Broglie period. The question is exactly how particle 
production compares with the electromagnetic loss rate. If the particle production rate 
is great enough, the decrease in particle number at some energy due to radiation can be 
completely compensated for. We estimate the magnitude of this effect for the specific 
form of fluctuating velocity above. In this very specific model we show that particle 
production is negligible. However, this is actually an artifact of the idealized choice of 
instantaneous switching between positive and negative fluctuations. We illustrate that 
in a more realistic model the amount of particle production is wildly dependent on the 
shape of the fluctuations. This makes it extremely difficult to derive any generic model 
independent constraints. 

4.1. Instantaneous random walk model 

We have assumed throughout this work a fluctuating velocity in the z-direction and so 
we need a model that gives such behavior. We will then show that tremendous amount 
of particle production can occur. The simplest example is that of a massive scalar fieldjl 
(p minimally coupled to a background metric of the form 

gaf3 = VaP + hafsit). (25) 

* We do not include more exotic scenarios like doubly special relativity. For an introduction to DSR 
see for example |22j . 

jj Our general conclusion, that particle production can be important, applies regardless of the bosonic 
or fermionic nature of the fields. Hence we will study a massive scalar field for simplicity. 
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hap{t) is a fluctuating term in the preferred frame tliat is dependent only on time. We 
assume that h^z is the only non-zero component as our choice of velocity fluctuations 
has been in the z direction. Furthermore we choose hzz{t) to be non- zero only in some 
finite region < t < T so that we have well defined asymptotic in and out states. The 
equation of motion for the scalar field in this background is 

^ --da (V^9''^d^<p) + = (26) 



which to lowest order in hzz is 

□0 + hzzit)dl^ - ^dthzzdt^ + m'ct) = 0. (27) 

The system is translationally invariant in z, so we can assume is of the form 
(j) = il){t)e^^^ . We then have the one dimensional equation for %l){t) 

+ El{l - hzz{t)) i) - dthzz dti) = (28) 

where Eq = p"^ + m? . 

The very specific case of the random walk scenario considered in section El can be 
implemented by choosing hzzif) to be 

hzz{t) = 5^ (t - ntf)e ((n + l)tj - t)A^ (29) 

n=0 

where = ±2fj. With this choice of hzz, dthzz = in the middle of the fluctuations 
(i.e. not at the points t = ntf). The damping term vanishes in these areas and (j28|) 
becomes a wave equation. To lowest order in the small quantities hzziTTt/p the group 
velocity Vg = dE/dp is 

.2 



v, = l-—,±Vf (30) 



m 
2^2 

which shows that this model realizes the random fluctuation scenario. The question of 
particle creation then becomes a matching problem between the boundaries of the n 
fluctuation regions as we transition from one wave solution to the next. 

The problem can be made much simpler by changing coordinates and eliminating 
the damping term entirely. We introduce a new time variable r] defined by 
f dt 

In ?7 time ()28p becomes 

i^{v) + E',(l-hl)4j{r^) = (32) 

where the double-dot denotes a double time derivative with respect to rj. In this 
form, ()32j) is simply a harmonic oscillator with a time dependent frequency. The key 
observation is that the frequency is quadratic in hzz- Since we are switching between 
ib2f J instantaneously the field equation in rj remains unchanged in the fluctuating region. 
Hence there is only particle production at the edges of the fluctuating region where the 



Constraints from cosmic rays on non- systematic Lorentz violation 14 

energy switches from E = Eq to Ef = Eq{1 — 3/4/i^^)^''^. Finally let us also observe 
from (p?T|l that for t < rj = t, while for t > T rj = t + r where r is a constant. The 
number operator in r) is therefore the same as that in t outside the fluctuating region. 
Likewise, the asymptotic in and out vacua are also equivalent, i.e. 

aJO >t= at\0 >r,= (33) 

where a^^t are the field annihilation operators. From these two observations it can be 
shown that the particle production calculated in rj is the same as would be calculated 
directly in t. 

We have established that in this model the only points relevant for particle 
production are the edges of the fluctuating region, where h^z instantaneously drops 
to zero. The amplitude for particle production in such a case can be calculated via the 
sudden approximation ^22j. At t = (when the particle enters the fluctuating region) 
the Bogoliubov coefficients are 



e- = -JJ^ + J^) (34) 




Eo 

With our particular form for Eq, Ef we have to lowest order in hzz 

a = 1 (35) 

/5 = ^hl- (36) 

When the particle exits aX rj = T + t this process happens again. The Bogoliubov 
coefficients have the same form as (jH^ but with, but with Ef and Eq interchanged and 
an additional phase e'^'^^'^~^'^\ The overall rate of particle production for an incoming N 
particle state can be calculated to be 

N,,, = NUl + P^htz) (37) 

where p is an 0(1) coefficient that depends on the phase. Since \hzz\ = ^Vf << 1 the 
amplitude for particle production in this current model is negligible. Hence the energy 
loss rate from radiation can actually be used to set a constraint on the size of Vf. 

4-2. Finite time model 

The lack of significant particle production in ()37p can be traced back to the fact that 
in the random walk scenario only the end points of the fluctuating region contribute 
to particle production. This in turn is a result of the instantaneous fluctuations. 
Realistically, of course, one might expect some finite time At] over which hzziv) changes. 
Over the course of the change there will be particle production. To show the model 
dependence, we return to (jH^ which applies even if the form of hzziv) is not that given 

by jir 
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Consider now one fluctuation, where hzz changes from 2vf to —2vf over time 
Arj. We define the deviation from 2vf by 2Ahzz = h^z — 2vf. We further assume 
Af] « Tjf where r// is the fluctuation time tf in rj coordinates. The hmit At] — is 
therefore the random walk scenario above, where every rjf time the fluctuations change 
instantaneously. If the amount of particle production from one fluctuation is small 
(which we show a posteriori) then we can calculate it via the background field method. 
We first rewrite (jH^ in terms of vj and Ah^z 

^ + E'^^l) - 3E^Ahzz{Ahzz + 2vf)tlj = (38) 

where E'"^ = Eq(1 — Sv'j). We define ^ = + ipi where i/jq = e"*'^''^ satisfies the 
equations of motion when Ahzz = 0. We are assuming there is only a small amount of 
particle production so ipi « ipo. To first order in the small quantities Ahzzy'ipi (jHH|) 
becomes 

-ipi + E'^V^i = 3E^Ahzz{Ahzz + 2vf)e-'^'''. (39) 

Equation (jH^ is now a harmonic oscillator with a time dependent force. The 
particle creation from such a system is well known (c.f. ^M,)- As an example, assume 
that Ahzz is linearff in rj, i.e. 

Ahzz = ^{v~Vo) (40) 

between times t]q and rjo + At] and zero elsewhere. Assuming ipi initially is zero, the 
outgoing number of ipi excitations after a single fluctuation is 

iV(^i) = 2E',Ar]^ (41) 

and the total expectation value of the particle number is l + N^il^i) since -00 was a single 
particle state. In the random walk scenario, tf was set by the de Broglie wavelength 
of the particle. If we set r]f by the same method, rjf = 1/E', and take Arj to be rjf/m 
where m is an integer then to lowest order N{ipi) = 2vj/m. 

For an individual fluctuation the amount of particle creation is actually small. 
However, the net effect of fluctuations is cumulative. Significant particle creation occurs 
when this cumulative effect becomes greater than or equal to the number of initial 
particles (which we chose to be one in the calculation above). After I fluctuations, the 
total number of particles created can be estimated by Ncum = {l+2vj/m)^ ^ l+2Ivj/m. 
For our previous constraint of = 10~^'^ and a proton of energy 10^^ eV travelling a Gpc 
this works out to be Ncum ~ 10^^/m. Even if m is of order 10^, which implies that the 
fluctuations switch over a timescale 1000 times shorter than the de Broglie wavelength of 
the UHECR protons, there would still be a magnification of 10^^ in the number of protons 
seen given an initial source flux. Such an increase in population is obviously unphysical, 
however our point is simply that particle creation effects can completely overwhelm any 
radiative losses and must be taken into account when deriving constraints. Since the 

ft This could be done by a more complicated functional form for hzz(t). Our point here is not to 
condone any particular model time dependence, but rather to show that particle production cannot be 
ignored. 
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amount of particle creation is model dependent little can be said about the expected 
population of UHECR protons in a Lorentz violating fluctuating background without a 
specific and concrete model for how quantum gravity might induce such fluctuations. 

5. Conclusion 

We have investigated in this paper a scenario with a fluctuating Lorentz violating 
dispersion relation from an effective field theory context. We have found two significant 
new effects. The first is that the low frequency radiation loss from such models is 
less than was calculated using threshold analyses. The actual constraint we get from 
UHECR protons is that the amplitude of the velocity fluctuations is bounded by 10~^'^. 
However, we also have shown that almost any constraint in a fluctuating background 
should be viewed with skepticism due to particle production from the time dependent 
background. For an initial one particle state, we find that for our constraint on Vf, 
there are of order 10^^ outgoing particles, completely overwhelming any threshold type 
analysis. The particle production is very model dependent, so it seems unlikely that 
any generic conclusions about fluctuating dispersion relations can be derived. It may 
be, however, that the particle production itself can be used to set limits for specific 
quantum gravity models. A similar approach has been used for causal sets in [7]. 

A limitation to our analysis is that we have strictly stuck to effective field theory. 
One could, of course, postulate that in some quantum gravity model such particle 
production does not happen. However, then we are well outside the realm of known 
physics and one must construct a new method for calculating reaction rates, etc. that 
yields no particle production for fluctuating backgrounds but also matches what we 
observe about relativity and the standard model. We know of no such framework, and 
hence we conclude that at the present time constraints from models that postulate 
fluctuating dispersion cannot be reliably extracted. 
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